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ABSTRACT 

In a finite-3 plasma ion cyclotron radius dispersion 
forms a trailing wave train for a perpendicular fast shock. 
Collisionless dissipation is provided by the three wave 
decay of the wave train into very oblique fast and parallel 
Alfven waves. . Particle thermalization results from Landau 
damping of oblique fast wave turbulbnce. The shock damping 
length to three wave decay is many ion cyclotron radii.' 
Undamped Alfven turbulence should persist far downstream 
from the shock. 
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1.0 INTRODUCTION ' 

+ ±2 

The high-3 (B = Suin' /B ) , truly collision-free plasma which 
permeates the interplanetary medium and most of the magnetosphere has 
stimulated considerable theoretical and experimental research on the 
structure of high-B collisionless shocks [Kennel and Sagdeev, 1967a, b; 
Tidman, 1967; Kennel and Pets chek, 1968; Fredricks et al. , 1970], 
Theoretical interest centers on determining the macroscopic or fluid 
behavior of the plasma flowing through a shock and specifying the 
collisionless plasma turbulence which provides microscopic dissipation. 
Two distinct approaches to this problem have evolved in the shock 
literature: a fully turbulent shock model in which collisionless 

dissipation dominates the structure [Fishman et al., 1960; Tidman, 1967; 
Kennel and Sagdeev, 1967a], and fluid shock models in which plasma 

•X. 

dynamics are governed by dispersive hydromagnetic equations and weak 
collisionless dissipation [Sagdeev, 1966]. The large parameter range 
and high variability of space plasmas permits ample room for application 
of bo^h methods. 

In this paper the fluid approach is employed to investigate a fast 
shock propagating strictly perpendicular to the ma^etic field in a 
high-B plasma. The shock structure is a nonlinear wave train with ion 
cyclotron radius (ICR) oscillation lengths which trails in the downstream 
flow. This nonlinear wave train decays by the three wave mode coupling 
process into oblique fast and parallel Alfven waves, thus relinquishing 
ordered shock oscillation energy to incoherent wave turbulence. The 
irreversible energy exchange between waves is the collisiohless 



dissipation needed to spatially damp ’.vave'train oscillations and complete 
the shock transition. . • ' 

Kennel and Sagdeev [1967b] first pointed out that finite- B pei^ien- 
dicular shocks should possess ion cyclotron radius wave train structure. 

As recognized by MacMahon [1968] and by Fredricks and Kennel [1968], 
however, Kennel and Sagdeev used an incomplete set of ICR fluid equations, 
and hence obtained the incorrect sign of ICR dispersion for the fast 
mode. (They found the ICR wave train leading in the upstream flow.) 
Consequently the three wave decay of the wave train into perpendicular 
magnetosonic turbulence proposed for the shock dissipation is invalid. 
However, the physical discussion of wave train formation, the three wave 
decay instability, as well as the basic analytical techniques emjJloyed 
to investigate wave trains and turbulent dissipation presented in Kennel 
and Sagdeev [1967b] are essentially coi*rect and are recommended to the 
interested reader. Therefore in this paper there is little need to dwell 
extensively on the physics of wave trains'or to expound on calculational 
details. 

Section 2.1 presents and briefly comments on the set of fluid 
equations, Chew-Goldberger-Low hydromagnetics with first order ion 
cyclotron radius corrections, used in the analysis. The linear dispersion 
relation for the perpendicular fast wave is derived. For. wavelengths 
comparable to the ICR, the wave propagates dispersively wfth phase speeds 
less than the long wavelength hydromagnetic fast speed. 

In the shock interior the fluid transition between upstream and 
downstream flow states is organized into a train of coherent or laminar 
oscillations which are weakly damped by collisionless dissipation. 

These, wave trains consist of dispersively propagating waves whose phase 
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speed matches the local fluid velocity; tlie relevant standing waves must 
be on the same branch of the dispersion .relation as the v/aves which form 
the shock. Since IGR dispersion slows tlie fast wave speed, the wave train 
with ICR oscillation lengtlis trails in the downstream shock flow. In 
Section 2.2 a differential equation for the flow velocity through the shock 
is derived. Here ICR dispersion balances the nonlinear hydromagnetic 
steepening of the shock and establishes a steady shock flow. The scale 
length of the leading velocity gradient in the shock and the dov/ristream 
oscillation wavelength is determined from the differential equation in 
Section 2.3. The ICR wave train is found to be a valid description of the 
shock structure only if the shock is weak, i.e., if the Mach number is much 
less than two. Hence the present analysis does not apply to the strong . 
bow shock interaction near the subsolar region of the magnetosphere. ICR 
fast shock str,ucture might occur in the weak shock region on the magnotor- 
sphere flanks or for shocks propagating in tiie interplanetary medium. 

Section 3.0 analyzes the three wave decay of the nonlinear perpendic- 
ular wave train into oblique fast and parallel Alfven waves. The matrix 
elements coupling the three waves arc derived in Section 3.1, following 
the standard techniques [see Galeev and Karpman, 1963; Kadomtsev, 1965; 
Sagdeev, 1966; Sagdeev and Galeev, 1969]. The growth rate of the decay 
instability is estimated in terms of the wave train oscillation amplitude 
in Section 3.3. The Landau or transit time damping rate of the oblique 
fast waves is shown to be less than the decay growth rate provided that B 
does not greatly exceed unity. In Section 3.4 the wave train spatial 
damping provided by three wave decay turbulence is found to be many ICR 
scale lengths. Since Alfven waves are not Landau damped, Alfven turbulence 
generated by the decay instability should persist far downstream from the 


shock. 
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2.0 ION CYCLOTRON R>\DIUS WAVE TRAIN FOR PERPENDICULAR FAST SHOCKS 

This section briefly reviews previous vvork on the ion cyclotron 

radius (ICR) perpendicular fast shock wave train [see Kins inger and Auer, 
1969; Goldberg, 1970], and makes specific the problem at hand. The time 
dependent Chew-Go Idberger-Low (CGL) hydroma gnetic equations correct to 
first order in the small ICR expansion are presented in Section 2.1; 
the perpendicular fast mode dispersion relation is determined. Following 
standard analytical techniques, a differential equation for the time 
independent wave train structure is derived in Section 2.2; the wave 
train properties are examined in Section 2.3. 

2.1 THE CGL-ICR FLUID. SET 

For time scales long compared to the ion cyclotron period and space 
scales longer than the ICR, plasma dynamics are often reasonably v;ell 
approximated by moments of the Vlasov equation. Of course, phenom.ena 
which depend on details of the particles’ distribution functions, such as 
Landau and cyclotron resonance, and heat flow, cannot be treated by fluid 
equations. By systematically expanding the Vlasov moments a set of fluid 
equations with first order ICR corrections can be obtained. to describe 
motion hoth- parallel and peipendicular to the m.agnetic field. The formula - 

■ ' V , ■ " ■ ■ 

tion employed here is due to MacMahon [1965], and the interested reader 
is referred to his paper for details. A coihmon difficulty with equations 
derived from a moment hierarchy is that they fail to close, i.e., lov/er 
order moments depend on higher order ones. This difficulty is avoided 
for the degenerate case of shock propagation perpendicular to the magnetic 
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field [MacMahon-,,. 1965] and for small amplitude linear v;aves. The consid- 
erations of this paper are restricted to tliese cases and, tlierefore, the 
fluid equations presented below are a truncation of the full set. 

The equation for the continuity of mass flow is 



V*(pv) = 0 


( 2 . 1 ) 


where p is the mass density and is the fluid velocity. The 
dynamical force equation is ’ . 




J X JB 

+ ^ 

c 


( 2 . 2 ) 


where is the pressure tensor summed for ions and electrons, and 

the superscript (1) denotes retention of first order ICR contributions. 
J is the current density, ^ the magnetic field, and c is the 
velocity of light. Gaussian units are used. The relation between J_ 
and B is specified by Ampere's equation ^ 


The electric field 


V X B = — J . (2.3) 

— c — 

E is governed by the particularly simple Ohm's law 


V X ^ 

E + — — =0 . (2.4) 


The neglect of contributions to Ohm's law from finite ion and electron . 

• • , A ■ , : 

inertia, and resistivity is a further restriction on the calculations 

± ±2 

performed here, hlien S = 8nP /B > 1 , inertial dispersion is a 
smaller effect than ICR dispersion. Faraday's law 

3B . , 

= - I ^ . ( 2 . 5 ) 

conpletes the Maxwell set. 
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For shock investigations an equation expressing the conservation of 
energy is required 


3 


' 2 r \ , . . 2 


• 

r 2 p^i) pf^^ 

pv" . ptn ♦ 1 pO) . f 

+ V • 

DV 


2 J. 2 »» 8ir 

V J 



[2 p 2p J 


* vpfl) * - . c,i«) * 


4ir 




= 0 


( 2 . 6 ) 


and are pressures perpendicular and parallel to the 

magnetic 'field. Subsidiary equations for and are needed 

below, and will be written here in a form appropriate only for perpendicular 
shock and linear wave propagation 


dP 


(1) 


~ (v-v + v^*v) + = 0 


(2.7) 


1 

2 dt 


( 1 ) 


+ P. 


CD 

U 


|vv 






= -0 


( 2 . 8 ) 


where is the perpendicular gradient operator. Contributions to 

(2.6), (2.7), and (2.8) from zero order heat flow along ^ have been 
neglected; these heat flo^v’ moments depend on the exact shape of the. 
distribution function and are not determined by the moment equations. 
Their neglect represents the closure of the moment hierarchy. 

The quantities q'*’.^^^ and are the first order ICR heat 

flows perpendicular to B . Only is needed, and subject to the 

4 b ■ 

same restrictions appropriate to (2.7) and (2.8), takes the form 


q^CD ^ 


2B X r 




A 


.p(D 


VP 


(2.9) 


t 
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is the zeroth order ion pressure; = eB/M^c is the ion g>"ro- 
frequency; e is the electronic charge and the ion mass.- In (2.9) 

terms involving heat fiov; and fourth order moments have been neglected; 
these terms can be significant for plasmas whose distribution functions 
deviate greatly from MaxAvellian [HacMahon, 1965, 1968]. 

Since the factor 1/fi^ appears in (2.9), is formally of 

first order; hence, formally, could be replaced by , 

whereupon it is easily shown that 3^^^^ vanishes. ICIi dispersion, 
however, couples different degrees of freedom, and the lowest order ICR 
corrections to long wavelength hydromagnetics are products of two ICR 


terms. Therefore in (2.9) V'P''*-' must be retained since 


Tnt.c-'i- •v'r.'f-ri'iT.ria ( i ) 

contributes to, and is actually the dominant contribution to, ICR 
dispersion [MacMahon, 1968]. 


If ^ is assumed to be in the z direction, the only components of 


,( 1 ) 


?'*■' required below are p^^^ and P^^^ ; these are 
■ ^ XX xy . 


CD ^ p(i) _ > % 


P 

XX 




2Q^ dx 


,( 0 ) + 


,(1) _ •A du 


xy 


2fi^ dx 


( 2 . 10 ) 


( 2 . 11 ) 


u and Vy are the fluid velocities in the x and y direction, respectively. 

Before passing to the shock wave train, a useful exercise is to 
determine the linear dispersion relation for the perpendicular fast wave. 
Linearizing about a stationary, uniform equilibrium of infinite spatial 
extent and assuming harmonic perturbations of the form 6u(xyz,t) = 

6u exp[ikx-ut] , (2.1) - (2.11) reduce to the dispersion relation 



-lo- 


co , . „ 2 


pCO)t 

3 • 

n - - 

^ 4 - 2 

PCn 


2 2 
^ J 

«2a " 


( 2 . 12 ) 


Cp is the hydromgnetic fast speed, ~ [C^^/4frp) + 
and is the ion cyclotron radius based on the 

zeroth order ion pressure. The small ICR expansion of the fluid 
equations requires that i -J shorter y-'ave lengths can only be 

treated correctly by the kinetic theory [Fredricks, 1968], 

As found by MacMahon [1968] and by Fredricks and Kennel [1968], ICR 
2 2 

dispersion, ,k R^- term in (2.12), decreases the phase velocity, below 
the hydromagnetic fast speed. In the context of shock wave trains, 
therefore, the dispersively propagating fast wave can phase stand, 
co/k = u = flow velocity, only in the doimstream shock flow. (The down- 
stream flow velocity must be less than Cp by the shock evolutionary 
conditions [Kantrowitz and Petschek, 1966],) Hence linear wave theory 
predicts that the 3 >1 perpendicular fast shock wave train trails 
behind the shock leading edge. 


2.2 DERIVATION OF THE SHOCK WAVE TRAIN : 

-When dissipative processes are weak, the collisionless shock 
transition from upstream to downstream stationary states takes the form 
of a nonlinear dispersively propagating wave whose phase velocity 
matches the local flow speed. This wave train is described by a 
differential equation which balances dispersion and weak dissipation 
against nonlinear shock steepening. Since the analytical techniques 
for obtaining the differential equation are standard [see Sagdeev, 1966; 
Cavaiiere and Englemann, 1967; Kennel and Sagdeev, 1967b], only the 
highlights are presented. A small viscosity is included in the analysis 
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to simulate weak collisionless dissipation, a specific model of vduch 
is the subject of Section- 3.0, 

Wave train solutions to (2.1) - (2.11) are sought which are time 
independent in the co-moving sliock frame. The shock normal is assumed 
to be in the x-direction, and is perpendicular to the magnetic field in. 
the z-direction; plasma quantities are assumed to be functions of x, only. 
Equations (2.1) - (2.6) now possess a first integral with respect to x; 
in regions v/here the flow is spatially uniform, these integrals are the 
f^iliar shock Rankine-Hugoniot relations whidi connect the upstreajn and' 
do\mstream states. For (2.1) - (2.6), however, only five independent 
relations between the six dynamical variables are obtained; hence an 
additional conservation law is needed. 

Goldberg [1970] has pointed out that for perpendicular propagation, 
the equation for , (2.8), provides the required conservation law^. 

In the shock frame (2.8) becomes 


~ (P' 
dx ^ \\ 


u + 23 * 1 ^^^ ♦ x) = 0 


J. 


(2.13) 


nfll ^ 

Since q *x vanishes for a spatially uniform flow, (2.13) yields 
a P^^i^Vp = = const. ; for perpendicular shocks the parallel 

temperature is unchanged. 

Since viscous dissipation is assumed small, second order terms in 
the coefficient of viscosity, y , are dropped in the derivation; in 
addition, terms of order |i(R_^/L) , where L is a characteristic shock 
scale length, can be neglected. With these approximations, (2.1) - (2.6) 
and (2.13) reduce to the following wave train differential equation: 



- 12 - 


1 

4 + dx 


rp(03+ 

~n dx 


yu . d Au 


Pl^l 


dx 


= u. 


Mp -1 

3 '.2 . 

■j (Au) + Au 2 — 


(2.14) 


m: 


u, denotes the upstream flov;- speed, and M =‘u /C is the magneto- 

■ 1 ^1 

sonic Madi Number; Au = (u-u^) /Uj^ . 

The first term in (2.14) represents ICR dispersion and will yield 
an oscillatory trailing v;ave train. Tne first derivative term repre- 
sents v/eak viscous dissipation. The right hand side of this second order 
equation can be interpreted as a nonlinear driving force; setting this 
term equal to zero yields the Rankine-Hugoniot relations for u . An 
analogy between (2.14) and weakly damped particle motion in an anharmonic 
potential well is presented in Sagdeev [1966] and in Kennel and Sagdcev 
[1967b] . The next section discusses an alternative method for obtaining 
the wave train properties. _ '■ 


2.3 STRUCTURE OF THE ICR WAVE TRAIN 

The wave train differential equation is nonlinear since both the 
coefficients of the differential operators and the driving force are 
nonlinear functionals of the velocity. Equation (2.14), however, does 
possess tv/o stationary points located at the upstream and downstream 
solutions of the Rankine-Hugoniot relations [zeros of the right hand side 
of (2.14)]. A general technique for determining the qualitative features 
of the nonlinear equation is to linearize about the stationary points 
and examine the stability of the linearized solutions. For wave trains 
• shock transition requires that the perturbations be unstable (stable) 



upstream (dovvnstream) [see Sagdeev/ ,1966; Cavaliere and Englemann, 1907], 
The linearized solutions determine v.'hether the v;aye train oscillations 
occur in the upstream or dovmstream flow and provide estimates of the 
shock scale lengths. 

After some manipulation, tlie linearized version of [2. 14) for the 
perturbed velocity 6u becomes 


3 p 2 d^6u jju d6u 

r* dx' \w)» dx 




6u 


0 


(2.15) 


where the coefficients of 6u are now to be evaluated about either the 

2 

upstream or dovmstream flow. If terms of orders p are neglected, 
solutions to (2.15) in the form 6u exp(Xx) are 


X = 


-pu 




L 

2 ~ R. 


3p(0). 

4 


1/2 


(2.16) 


Since . Mp -1 > 0 upstream, (2.16) yields u decreasing exponentially 


from Uj with a scale height given by 


L 'V' 


3 

4 


p(0)+il/2 




1/2 


(2.17) 


Hence ICR dispersion determines the thickness of the leading edge. 

2 

Downstream M_ -1 < 0 , and solutions of (2,16) are damped oscillations 
r 

with wavelengths given by (2.17) evaluated about the downstream flow 
and damping length j . Tnerefore as predicted by the 

linear wave dispersion relation [Macraahon, 1968; Fredricks and Kennel, 
1968] ICR dispersion forms a trailing shock wave train. 



The first .order TCR-CGL fluid equations provide an accurate 

description of tlie shock structure if, and only if, the shock .thickness 

L « . Therefore from (2.17) the upstream Mach number is limited to 

values Mp -1 << < 1 , so that ICR dispersive shocks must be 

weak. Coroniti [1970] has shown that stronger high-B fast shocks 

steepen until thicknesses of order the electron* intertia , length C/o) 

2 2 

(w^ = 47rNe /M ) , are reached; for solar wind plasmas, « R^ . 

The shock structure is now described by a trailing length wave 

train. 

In 'summary, the ICR-CGL hydromagnetic equations describe plasma 

behavior over scale lengths long compared to R^. A closed set of 

Rankine-Kugoniot relations for perpendicular shocks exists with Tj^ 

constant across the shock. iVhen ICR dispersion is balanced against 

nonlinear steepening j a trailing wave train with scale lengths 

resolves the shock structure. This shock solution is valid only for 
2 

weak shocks, Mp-1 « 1 . In the above analysis weak collisionless 
dissipation was assumed to provide an irreversible shock transition. 

The next section considers a specific model, the nonlinear three wave 
decay process, for this dissipation. 



3.0 THREE WAVE DECAY. INSTABILITY OF A NONLINEAR PERPENDICULAR FAST V/AVCE 

3.1 INTRODUCTION 

A nonlinear wave, such as a wave train, docs not always constitute 
an equilibrium state since the oscillation amplitude is far from the 
thermal level. Two perturbation vvaves, either thermal or nonthermal, 
fluctuations, can couple to or scatter off the nonlinear wave spatial 
gradients [Galeev and Karpman, 1963; Kadomtsev, 1965]. If these mode 
couplings resonantly extract energy from the nonlinear wave, thus 
amplifying the perturbation waves, the nonlinear wave amplitude can 
decay faster than the rate set by particle-particle collisions. A 
nonlinear wave train, therefore, can parametrically excite incoherent 
wave turbulence which, in general, will not phase stand in the shock 
flow, and which provides a form of irreversible collisionless dissipation 
that damps the wave train amplitude. Particle distributions adjust by 
damping the turbulent waves, eventually relaxing to a uniform downstream 
state. 

The three v/ave decay instability is possible only if the wave 
frequencies and wave vectors of the interacting modes satisfy the following 
resonance conditions [Sagdeev and Galeev, 1969] 

Uj = Uq + 0)2 (3.1) 

k 1 = + k 2 (3.2) 

0 )q, refer to the nonlinear wave, and 0 )^^, and 0 ) 2 , ^2 
two perturbing waves.- In analogy with quantum mechanics, (3.1) and 



(3.2) express the conservation of wave energy and wave mbnientum- in the 
/ ■ 

interaction. In addition to (3.1) and (3.2) the matrix elements coupling 
the three weaves must be non-zero and have the appropriate sign for 
instability. 

A general theorem governing possible decays is that three waves on 
the same branch of the dispersion relation can mode couple only if to 
increases with increasing k (see Sagdeev and Galeev [1969] or Kennel 
and Sagdeev [1967b] for proof) . Therefore, from (2.12), decays betv.-een 
three perpendicular high-3 fast waves are disallov/ed [MacMahon, 1968]. 

The fast mode, however, can unstably interact with waves on other branches 
of the dispersion relation provided that the decay is not forbidden by 
polarization restrictions.' 

This section considers the decay of the perpendicularv ICR wave train 
into slightly oblique fast waves and Alfven waves propagating parallel 
to the magnetic field. The motivation for treating this-particular case, 
is that Alfven waves are undamped in high-3 plasmas [Barnes, 1967], and 
therefore could possibly be observed in the downstream flow by space- 
craft. Other decay modes are not considered but could contribute to the. 
total shock dissipation. The matrix elements for the interaction are 
derived in Section 3.2. Tlie growth rate for the decay instability is 
determined in Section 3.3. Section 3.4 estimates the wave train damping 
length. ■ 

3.2 THE MODE COUPLING MATRIX ELEMENTS 

2 

Since the ICR shock must be weak, Mp - 1 « 1 , the do^^nrlStream 
wave train oscillation amplitude is small. Therefore locally the wave 
train should approximately obey the linear wave dispersion relation, 


s 
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and the wave polarizations canbe interrelated hy linear theory. For 
the perpendicular ICR shock wave train, the polarizations are • 6V , 6V , 

^ y 

6p, , and ; the temporal and spatial dependence is of harmonic 

form exp [ikj^ x-iw^t] . The slightly oblique fast mode is denoted by 

° Cl) 

Xi» ki* ^1 ’ oscillates as exp [i (kj^ x + kjj z) - ico^t] . 

Similarly the Alfven wave polarizations are V , b , which vary as 

^2 ^2 

exp[ik.. 2 - ito t] , 

2 ^ 

From C3.2), the wave vectors of the three interacting modes must 

satisfy k. = k, E k, and k„ = k,, = k,, . , For slightly oblique 

0 ^ "2 

propagation of wave (1), k^ » k,j , an assumption whidi permits 
considerable simplification in determining the matrix elements coupling 
the three waves. 

The calculation follows the method of Sagdeev and Galeev [1969]. 

The dispersion relation for each mode is determined from Eqs. (2.1) - . 
(2.9) retaining the resonant nonlinear coupling of the perturbation waves 
to the large amplitude wave; the resonant couplings are of the form 
6V V for the Alfven wave and 5V* V for the oblique fast mode 


X X 




(* denotes complex conjugate). Other nonlinear terms such as V or 

■■■ ^1 

V V are neglected. The linear wave polarizations are then used to 

*1/2 

eliminate all the variables in terms of V , V , and 6V . The wave 

^1 ^2 . 

polarizations are nonlinearly coupled by matrix elements, the relative 
signs of which determine the stability of the interaction. Since the 
analytical procedures are standard, the calculation is outlined only for 
the Alfven wave; the result for the fast mode is simply stated. Further- 
more assuming k » k , ICR terms of order k R « 1 contribute 
negligibly compared to k^_^R_^ < 1 terms and are dropped. 



f 


Substitution of the above polarizations into the y-component of 
Ohm's law, (2.4) and (2.5), and use of (3.1) and (3.2) yields the 
following equation for b 


exp[ik|jZ-U2t] 




C^V*b^^+Vy 6B*) exp[ikyZ-iw2t] 


+ nonresonant terms 


is the equilibrium magnetic field strength. The nonresonant terms 
vanish after the phase averaging. The lowest order linear polarizations 
for ^^ 2 ’ ^y are given by 


6B = 
z 


A X 0 


6V = i 

y 


k.^R 6V 

i + + X 


V B- 

i Xj 0 


k.^R V 

Ji + + X, 

i ^ 2 2 — - 

2[Wy-(k/G^/u)^)] 


2 1/2 

where = (B^ /4Trp) , the Alfven speed. Substitution of (3.4) 


(3.7) into (3.3) yields 


. k V B^ k k.Sl B-k ,^R ^6V*V 
H X2 0 . IJ A + 0 4 + . X Xj 

^^2 " ^ ^ ^ 20 )^ 10 ^ 0)2 


k 2c 2 

o)j2[l-(kjj^C^2/o)^2)^ 
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• 2 2 2 

Note that since' the oblique fast mode has w, h.C^ . k„ C /w, . 

1 F ’ A / 1 

2 2 

kjj /k^ . « 1 , and the nonlinear contribution to (3.8) is small. 


Performing an identical analysis on the y momentum equation yields 


for 'V 




ik B b 
II 0 


4irp 


+ 1 


0 


k B-b 
11 0 Xl 1 

' 6B* 
z 

' 

6p*‘ 

o 

Q. 

■ 

Po 

J 


+ ik. (V 5V* - V 6V*) 
X Xj y Xj x" 


(3.9) 


The last two nonlinear terms, arise from- V*VV . The perpendicular fast 

wave polarizations satisfy 5(B^/p) = 0 ; hence the nonlinear term 

proportional to b vanishes identically. Substituting (3.5), (3.7), 

2 2 2 

and (3.8) into (3.9) and neglecting terms of order k^^ /oi^ « 1 , 
yields the following simple result 


w. 


k 2 ^ 
^11 


2 2 
k.fi k, R ^ 
+ 4 + 


L 

W, 




6V*V 

X X, 


(3.10) 


The coefficient of V is just the linear Alfven wave dispersion 

y 2 ' 

*■ 2 2 
relation. Note that the nonlinear terms are proportional to • 

A more convenient and transjjarent form of (3.10) is obtained by 

reformulating the above analysis in a Hamiltonian representation [Sagdeey 

and Galeev, 1969]. The time dependence of the wave amplitude is now 

assumed to have the form V = V (t) exp[.ik z-iu„t] where V (t) is 

^2 ^2 r ^ '2 

considered as slowly varying. Use of the linear dispersion relation 

for then permits (3.10) to be written as 


3V ■ (t) 


8t 


2 2 

kj^n.k.^R 

2 03 , 


K. K p 03- - 

-- -- 2 - — 03 - 03 , 6V*V 

Iq03j^ 052 L J ^ ^ X Xj^ 


(3.11) 
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v.'here ,^'^q = ^1~^2 used. « . . ; 

An identical analysis for the oblique high-B fast mode yields for 

V (t) 

1 


_1 li k X + + 


6V V 


(3.12) 


2 2 2 

where terms of. order k,, C /w « 1 and k„R « 1 have been 

•I A i 13 + 

neglected. A similar equation for the time rate of change of 6V 


could be obtained. 


3.3 DECAY INSTABILITY GROWTH RATE 

Now consider. the simple case where the wave train amplitude greatly 

exceeds both perturbation an^jlitudes, 6V » V ,V . If the growth 

X Xj y2 

time for the instability is long compared to a flow time across a wave 

train scale length, the time dependence of 6V^ in (3.1!) and (3.12) 

can be neglected in the determination' of the initial decay instability 

growth rate. Of course after many e-folding times of the instability, 

6V^ will be reduced to the perturbation amplitude level; however, by 

this time the shock transition will be completed. Note that the nonlinear 

limit cycle of the decay instability, i.e., when all three modes share 

equal amplitudes [Sagdeev and Galeev, 1969] , need not be considered since 

the perturbed waves are convected downstream, out of the shock I'egion. 

If V (t) and V (t) vary as exp[yt] , (3.11) and (3.12) 

^1 ^2 

reduce to • 


2 

Y = 


k \ \ Ik ^ 


4|a)Qa)jW2 


[^ -^] “l“2 


2 . 


(3.13) 
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Instability results only if sign <0 , or from (3.1) if- 

Wp > + 1^2 1 . .Hence the- initial large amplitude wave decays to. 

waves of lower frequency or, by the quantum analog)’’, the decay occurs 
only to lov/er energ)’’ states. 

To further determine the growth rate, hj, , , and must 

be estimated. Setting co^ = ^ ^ ' *^0 ” C^*!) 

and (2.12) w, becom.es w, = M„ ^ C., - k.r. = k, C„ [1 - (3/8o) (Fp^‘"/C^^)k /'R ^ 
i i. r.., r "A r r « + 

2 2 ^ 

where k, R <1 was assumed. .Note tliat w, /k. < u« ; hence the 

perturbed .waves are converted downstream. From the Ra'nkine-Hugoniot 

relations for a perpendicular sliock [Anderson, 1963], M_ is given by 

. 2 


M, 


r (M: -1) + 3 ql/2 
1 

4(mJ, -1) + 5 


Mr 


(3.14) 


if Mp -1 « 1 . Eliminating H„ then yields 
1 ^2 


2k C pC0) + 

.. ,„2 , , . 3 A 2„ 2 ^ A v-j 

-1^- -(M -1) + 4 — T2 \ 

A F 1 PpCp 


p(0) + 

3 X , 2„ 2 

4 “TT 


(|2gj2/w^|) 


(3.15) 


0 ) 2/031 


The wave 


as an order of magnitude estimate for k jj and 

2 

train amplitude is estimated as 6V^ ')/ (dp/p^) (oj^/kj^) Cp(Mp -1) 

0 0-' 21 
Cp(Hp -1) . Substitution for o)q, <5V^, and 1 20 ) 2 / 031 1 into (3.13) 

yields as a very rough estimate 
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Y a. 


(0)-v 


4 p. C. 


(5.16) 


0 r 


Tlie above calculation assurasd that oblique fast aiid Alfven v;aves 
were present in the noise ..fluctuation spectiaim and could be parametrically 
amplified, by three wave decay. Oblique, fast waves ,.. but. not Alfven v/aves, 
have resonant Landau or transit time damping interactions v.dth thermal 
ions and electrons [Stix, 1962; Barnes, 1967]. If the nonlinear growth 
rate (3.16) is less than the Landau damping decrement, then little or no 
energy Can be coupled into fast waves by nonlinear mode coupling, and the 
three wave decay instability probably does not provide the shock dissipa- 
tion. The oblique fast wave Landau damping rate is approximately [Stix, 
1962] 




(3.17) 


2 ± ~ 

where C.^. = 2T , the specie thermal speed. Substitution of (3.15) 

for k , use of the dowT.stream phase standing condition 0 )^ = = 

~ T estimate u and 

2 2 

k. R in terms of M , and elimination of M„ by (3.14), yields for 
i + ^2 ^2 

Iy,"| 


. + C. k,C„ 

I Vi = 

^ (M^ -1) 


exp 


4C 

^^A 

Tvrr 


L c, (Hp-D-J 


(3.18) 


The ratio of the nonlinear growth to the linear damping rate then becomes 


,(0) + 


C. 


Iy, 




2 2 2 2 
kj^^R^^(Mp^-l)'' exp 


AC. 


2 2 2 

c/(Mp -1) 

- ^1 


(3.19) 



Since C v4'I^'/M l/</^ » 1 unless; 3 M^/M' > , electron 

Landau damping does not suppi'ess the nonlinear decay instability.. Ion 
Landau damping is also small provided the exponential factor is large, or 

Mj. -1 « 2//fi . Hence three v.-ave decay is an effective dissipation 

1 

process for the perpendicular fast shock if 3 1 , but becoBies 

. ■■ ■ ‘ i ^ 

considerably weaker for higher 3 • 


3.4 ESTIMATE OF THE WAVE TRAIN DAMPING LENGTH 


The three v/ave decay instability extracts enei-gy from the ordered 
oscillations of the v/ave train and generates incoherent wave turbulence 
in the downstream flow. Hence the v;ave train should pei'sist for a length 
L„ 'w M Cp/y.. before being damped to the fluctuation level. Using (3.16) 
Ljj . can be estimated as 


, 4 foS' 

T T V3 pW* 


3 3 

T -v. 


(3.20) 


From the dovnistream phase standing condition, ~ Mp and (2.12) 

and (3.14), is estimated as ~ y ’ 


Lp then becomes of order 


L 
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f p(0)r 


Ipo"" 


1/2 


R. 


F J 


Ol -1)^'" 
b 


(3.21) 


For Mp.-l « 1 , the wave train should exist many ion gyro-radii in 
the downstream flow. ' 

If no other particle dissipation processes, such as current 
instabilities, are involved in the shock structure, shock heating of the 
particle distributions occurs primarily by. Landau damping of oblique 
fast wavj:s. The particle thermalization length should be of order 


I 
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, or substituting- (3,18). 


L - vi 

T 2 




7rpC, 


1/2 


„1_ -1)^/^ R exp 

0-Ca 4 


4C 
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LC^^(Mp^-l)^ 


(3.22) 


•T* 

Note that by (3.19) L.^ > ; hence particle dissipation persists 

further downstream than does wave train dissipation. 

i 

Since Alfven waves are undamped, after many or L.j,~ lengths 

downstream the wave spcctnim w-ill be dominated by transverse magnetic 
oscillations. 


3.5 SUM>-1ARY 

The parametric amplification of oblique fast and- parallel Alfven 

v;aves at the -expehse bf the Wave train amplitude produces the shock - 

dissipation required by the Rankine-Hiigoniot relations. In the downstream 

flow coherent v/ave train oscillations gradually become disordered and 

finally damp to the fluctuation level. Oblique fast wave turbulence is 

Landau damped by ions and electrons, thus producing downstream. particle 

thermal! zati on., until only Alfven turbulence remains. The relaxation 

lengths for the v;ave train and particles are many tim,es R^ ; hence the 

complete shock transition . region .should be much broader than the initial 

gradient of the magnetic field at the leading edge. Landau damping of 
■ . ■ ■ 
oblique fast waves probably restricts the particular decay instability 

i ■ ■ 

dissipation considered, here to shock flows in v;hich 3 does, not greatly 
exceed unity. .The spatial structure of the shock is sketched in Fig. 1. 



s 


-25- 


4 . 0 DISCUSSIO N 

The three wave, decay model of .collisionless dissipation probably 
greatly oversimplifies the turbulent relaxation of tlie downstream shock . . 
flov;. Clearly; decays into other modes of hydromagnetic turbulence are 
possible and v.’ould also contribute to wave train dissipation. Hence 
the dariiping length determined from a single decay instability only . 
represents an upper limit to tlie shock length. 

Heat transport in tlie downstream flow might also be quite complicated. 
If particle thermalization proceeds by the linear Landau damping of 
oblique fast wave turbulence, the quasi- linear theory predicts that the 
resonant diffusion is entirely in parallel energy [Kennel and Englemann, 
1966],. The nankine-Hugoniot relations, hov.-ever, require that T^ remain 
unchanged across the shock; hence additional turbulent dissipation is 
needed to convert acquired parallel energy into perpendicular energy. 

For Tjj > Tj^ and B large, parallel propagating fast waves are 
unstable to the nonresonant fire-hose and to a resonant ion instability 
[Kennel and Sagdeev, 1967a; Kennel and Petschek, 1966]. Kennel and 
Scarf [1968] have shown that the ion resonant grow'th rates are exponentially 
enhanced if the electrons are also Tj~ > Xl anisotropic; the instability 

now can affect the .main bulk of the' ion distribution, thus greatly 

\ ■ , . 

increasing the turbulent dissipation. The downstream wave spectrum, 

therefore, is likely. to exhibit a variety of turbulent modes each attempting 
to drive the plasma toward a stationary configuration. 

The decay instability considered here predicts that Alfven wave 
turbulence should persist far downstream from the shock. Detection of 
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tliese waves , however, does not itself-.constitute a verificat-ion. of ..the '.. . 

proposed dissipation model since the fully turbulent shock theories of 

Kennel and Sagdeev [1967a] and of Camac et al. [1962] also predict 

downstream Alfven turbulence. Direct observations of the internal shock 

structure are required to differentiate betv/een the various tlieories. 

Typical solar v/ihd flows have moderate to higli Mach numbers 

Kip 4 to 10. Therefore in the spatial region surrounding the bow of . 

the magnetosphere where the shock is strong, trailing ICR wave trains 

of the t>’pe discussed here will not generally be part of the shock 

structure. Fredricks etal. [1970] have . reported that near the sub-solar 

region the bow shock is characterized by electron inertia lengths , C/w . 

+ 

which ax'e much less than R when 3 1 . Far out on the flanks of 

+ 

the magnetosphere, however, the shock strength is weak. Also shocks 
propagating in the interplanetary medium might be weak,. Hence the 
trailing ICR fast shock structure could be detected in these regions by 
high teleractr}'' rate satellites, . ._ 
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FIGURE CAPTION 

Fig. 1. Tne spatial structure of the magnetic field in the ICR fast 

shock is sketched against distance through the shock. Upstream 
is to the left, downstream to the riglit. At. tl\e leading edge 
B increases exponentially with a scale height L "w .■ V/ave 
train oscillations , also witli R scale lengths, trail in the 
downstream flow. The wave train decays into oblique fast. and 
-parallel Alfven waves by three wave mode couplings, thus 
providing irreversible collisionless dissipation, have train 
oscillations are damped over a length . Ion and 

electron thermal! zation proceeds by the Landau damping of 
oblique fast waves. Undamped Alfven waves should be observed 
far dov;nstream. 




